Effects of a Fluctuating Interface between a Superfluid and a Polarized Fermi Gas 
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Motivated by recent experiments in trapped Fermi gas with spin population imbalance, we discuss 
the effects of the quantum and thermal fluctuations of the interface between a fully paired super- 
fluid core and a fully polarized Fermi gas. We demonstrate that even if there is no true partially 
polarized thermodynamic phase in bulk, the interface fluctuation can give rise to a partially polar- 
ized transition regime in trap. Our theory yields a definite prediction for the functional forms of 
the spatial profile of spin polarization and pairing gap, and we show that the spin-resolved density 
profiles measured by both the MIT and Rice groups obey this function form. We also show that 
sufficient large fluctuation will lead to a visibly unequal density even at the center of the cloud. We 
hope this picture can shed lights on the controversial discrepancies in recent experiments. 
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Phase separation is an ubiquitous natural phe- 
nomenon. Its root is first order phase transition. For 
example, as a function of temperature (T) and mag- 
netic field (if), an easy-axis ferromagnet phase sepa- 
rates into two oppositely magnetized phases when (T, H) 
are fined tuned to the boundary of the first-order tran- 
sition (H = 0, T = T c ). On the other hand, as long as 
T < T c phase separation occurs without needing fine tun- 
ing when the magnetization, rather than magnetic field, 
is fixed. In the presence of phase separation the dynamics 
of the interface determines the low energy/temperature 
physical properties of the system. Recently experiments 
of ultracold Fermi gases near Feshbach resonance once 
again reveal phase separation phenomenon [l|, 0, El- I n 
these experiments, an external-imposed spin population 
imbalance frustrates the strong pairing interaction due to 
Feshbach resonance. As a result, the system phase sepa- 
rates into a core of superfluid with equal spin population 
and a shell of excess unpaired fermions. 

Although similar phenomena were observed, there are 
important difference between the findings of two major 
experimental groups at MIT and Rice University. Let 
Ni,i represent the number of spin up/down atoms, and 
P = (JVf — Ni)/(Nf + Ni) be the population imbalance. 
Here are the highlight of the differences. The MIT group 
finds (i) for small P, the spin-resolved density profile ex- 
hibits three regions: a superfluid core with n-\ = tij, a 
partially polarized (PP) shell with nj > rij > and a 
fully polarized region with — 0. The width of the PP 
shell is as wide as the other two regions (see Fig. 5 (a) of 
Ref.Q). (ii) The interfaces between different regions al- 
ways follow the equal potential contour as predicted by 
theories using local density approximation, (iii) As P 
exceeds a critical value P c ~ 70%, the equal density core 
disappears and also the condensation fraction vanishes, 
which is referred to as the Chandrasckhar-Clogston (CC) 
limit of superfluidity. In contrast, the Rice group finds 
(i) the width of the PP region is much narrower than the 
width of other two regions (see Fig.3(b) ofRef.(|)). (ii) 



At large P the interface does not follow the equal poten- 
tial contour, (iii) The equal density core exists up to the 
highest P studied (~ 95%), and there is no sign of CC 
limit. 0] One can expect the explanation of these discrep- 
ancies to shed light on a number of fundamental issues 
of this problem. For examples, what is the nature of the 
PP region, does it reveal a true quantum phase of the 
system, and what controls its width; does the CC limit 
for superfluidity exists, and if so, what is the mechanism. 

These experiments have stimulated much theoretical 
discussions in the literature. Most of the existing the- 
ories assume the system to follow the confining poten- 
tial adiabatically, so that locally it corresponds to a bulk 
phase with the chemical potential /i^j = /lot.i — ^( r )j 
where /^ot,i is a spin-dependent constant and V(r) is the 
trapping potential. The issue at hand is whether the PP 
state is a stable bulk phase near Feshbach resonance. By 
comparing mean-field energies, many authors predicted a 
narrow window of PP phase in the bulk phase diagram^ . 
However due to the approximate nature of the treatment, 
and the fact that even in mean-field theory the energy 
difference between different phases is tiny, it is hard to 




FIG. 1: Schematic of the interface flucuation. The dashing 
line is the static location of the interface which separates a 
fully paired superfluid phase and a fully polarized fermi gas. 
The solid line sketches an example of surface mode which 
distorts the interface. 
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know whether such a PP bulk phase really exists or not. 

In the following we make a bold assumption: the PP re- 
gion found in experiments is not a true phase, rather it is 
due to the fluctuations of the interface between the equal- 
population supcrfluid and the fully polarized phase. Our 
assumption implies that, for interaction strength rele- 
vant to the experimental systems, there is a direct zero- 
temperature first order transition between the superfluid 
and the fully polarized phases as a function of the Zceman 
field. Our proposal is motivated by two recent develop- 
ments, (i) Recently the MIT group found the presence 
of a full pairing gap in the PP state at P > P c where 
the superfluidity is destroyed 4]. This contradicts theo- 
ries which attribute the CC limit to depairing|5j. (ii) The 
importance of interface physics is highlighted by recent 
theoretical works which attribute the deformation of the 
interface from the equalpotential surface to the existence 
of finite surface tension]?], Q . Of course, interface physics 
involves not only statics but also dynamics . 

In the rest of the paper we analyze the effect of har- 
monic interface fluctuations between the superfluid and 
the fully polarized phase. In FigfT] we illustrate a snap 
shot of an allowed interface deformation where the en- 
closed volume is kept fixed. Surface fluctuations have 
been studied in early days of atomic BEC and fermion 
superfluid with equal population 9] . However since the 
system is consisted of only one single phase, such fluc- 
tuation is less important as it only affects the low den- 
sity region at the edge of the cloud. In comparison the 
effects of interface fluctuation are far more pronounced 
here which lead to a very precise and universal profile of 
spin polarization and superfluid order parameter across 
the two phases. 

Here are how we describe the interface, i) Motivated 
by the experimental facts 0, S|, we assume the equilib- 
rium interface S is a closed two-dimensional surface with 
cylindrical symmetry: 

S = {{x, y, z):x 2 +y 2 + e 2 (z)z 2 = a 2 }. (I) 

Here e(z) is an even function, its value is unity for z = 0, 
and decreases monotonically to o/zq < 1 as z approaches 
±zo, the north and south poles of S. In the following 
we shall use z and the azimuthal angle <fi as the co- 
ordinates for S, so that 



following differential equation 



x — ^/ a? — e 2 {z)z 2 cos (f>, y = 
a 2 — e 2 (z)z 2 sin0, z — z. It should be apparent that 
this is a locally orthogonal coordination, i.e., ds 2 = 
g ( j > {z)d4> 2 + g z (z)dz 2 . (ii) We describe the fluctuation 
of the interface from its equilibrium shape by a har- 
monic quantum elastic model. Let u(z, <fi, t) be the time- 
dependent normal displacement of the interface from its 
equilibrium position, the Lagrangian of our model reads 

V{z) 2 ^ (d,u) 2 i 



C= d4>dz y /g < t,g 2 
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,(2) 



where m is an introduced phenomenological parameter. 
The normal modes of the interface fluctuation obey the 




- V(z)u a = -~mu; 2 u a ,(3) 



where uj a is the normal mode frequency. In later analysis 
we shall expand u(z, (f>, t) in terms of the normal modes : 



i(z, <j>,t)=y^ j A a (t)u a (z, t 



(4) 



Note that all quantities in the above equation are real. 
In Eq. (|4]) the sum excludes the "breathing" mode, the 
mode whose associated distortion does not conserve the 
enclosed volume. Interestingly, despite the presence of 
g$(z) and g z (z) in Eq. ([2]) and Eq. ([3]), the functional 
form of the interface profile can be determined exactly 
along special directions such as the z-axisQ and any ra- 
dial direction in the z = plane 0]. 

For example, as the main result of this paper, we ob- 
tain the following function form for the spin polarization 
M = n.f — 7ij_ and the pairing gap profiles along a high 
symmetry radial direction (t): 



M(r,t) = - 



A(r,t) 



I + Erf 
I -Erf 



r - r (t) 

£(*) 
r - r (i) 



n p (r,i) 
AsM). (5) 



Here Erf(x) is the error function, n p (r,t) is the density 
profile of a fully polarized Fermi gas, and A B (r, t) is the 
pairing gap profile of the fully-paired superfluid, along 
the same direction in the same trap. In particular we 
use II ill 



n p (r,i) =n [ 1- 
A s (r,t) = A [ 1 



(r/R(f)) 2 f' 2 
-{r/R{t)f\. 



(6) 
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FIG. 2: Plots of Eq. © and Eq. © for different choices of 
i/R. We have use r /R = 0.3 for all curves. £/R = 0.03 
for green line, (,/R = O.f for blue line and £/R — 0.2 for red 
line. The black dashing lines are the profiles without interface 
fluctuation, which exhibit discontinuity. 
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FIG. 3: Fit the measured spin polarization profile to Eq. ([5]) 
and Eq. (|6]). The green line in (A) and the black line in (B) are 
the data of the Rice0] and MIT0] group, respectively. The 
black line in (A) and the purple line in (B) are our best fits. 
The fitting parameters for (A) are ro = 240/im, £ = 80/im, 
n = 2.97 x 10 n cm~ 3 and R = 800^m. The fitting parame- 
ters for the left and right panels of (B) are ro = 27/j.m, £ = 



23am, n = 2.058 x 10 i/ cm" J R 
£ = 35/xra, n = 2.645 x 10 12 cm -3 



= 70/im, and ro = 16/ira, 
R — 73[im, respectively. 



In Fig. ([3]) we fit the polarization profile measured 
by both the MIT and Rice group to the first line of 
Eq. ([5]) and Eq. ©. These fits are excellent! Recently a 
new scheme of tomographic rf spectroscopy invented in 
MIt[Ii[ gives the promise of measuring the spatial pro- 
file of pairing gap A(r), which can be used to test the 
second line of Eq.©- For the same sample, if one mea- 
sures both 8n(r) and A(r), and provided that both 5n(r) 
and A(r) can be well fit by Eq. one can extract the 
values of ro and £ by fitting these two spatial profiles 
to these two functions. A falsifiable predication of our 
phenomenological theory is that the values of r and £ 
extracted from these two spatial profiles, 5n(r) and A(r), 
shall be the same. This is because in both case it is the 
same interface fluctuation that gives rise to the smooth 
spatial profiles. 

Now we fill in the details of how to obtain Eq. ([5]) from 
the model defined in Eq. © and Eq. ©. We will first 
consider the effects of quantum fluctuation at T = 0. 
Later we shall generalize the result to finite temperature. 

Substituting Eq. Q into Eq. © and use the orthonor- 
mal relation between different normal modes 



^g4>g z d(i)dzu a {z, <j>)up(z, cj>) = S a ./3 



we obtain 



C 



\2 _ 2 A2 



(7) 



(8) 



Eq. ([5]) describes a set of independent harmonic oscilla- 
tors. The corresponding Hamiltonian is given by 



In writing down the first line of Eq. ([6]) we have assumed 
the temperature is sufficiently low and the trapping po- 
tential is harmonic. For the second line of Eq. © we 
have taken into account the universality of the pairing 
gap near Feshbach resonance (Iol|. 

An appealing feature of Eq. ([5]) is that the microscopic 
information about the interface (and temperature) en- 
ters through a single parameter £(t). In Fig. ([2]) we plot 
Eq. ([5]) and Eq. ^ for typical values of parameters. 
For ro ^> f, the density difference at the center of the 
cloud is exponentially small which reveals an equal den- 
sity core. For small £ (green line), one can have rather 
narrow intermediate partially polarized regime (like the 
Rice group's result). Wider partially polarized regime 
(like the MIT group's result) can be obtained when £ is 
given a larger value (blue line). When £ increases toward 
ro, hence the interface fluctuation becomes more severe, 
even the center of the core can become visibly partially 
polarized. Meanwhile the pairing gap Ao remains finite 
similar to the result of Ref. 0] • Clearly in this theory the 
partially polarized region is caused by the interface fluc- 
tuation rather than a PP normal phase with a vanishing 
(or smaller) pairing gap. 
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where II Q and Ap obey the following commutation rela- 
tion [A^ILJ = ihS a p. The ground state wavefunction of 
Eq. ([5]) is given by 



, (io) 



where Y\ excludes the breathing mode from the product. 
Here we explain why some directions (t in Eq. ([5]) and 
Eq. ([6])) are special. For a generic point on the surface 
specified in Eq. ([I]), the direction of surface normal does 
not agree with the vector connecting the center to the 
point in question. There are special points on the surface, 
for which these two directions agree. For the surface 
described by Eq. ([1]) these special directions are ±z and 
f around the equator at z = 0. Let (zo,<po) be such a 
special point, 



u (zo,4>o) = A a u a (z ,<t)o). 



(11) 
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Let us first calculate the probability density for the in- 
terface to move by amount n 

P{n) = fX^dA^irt-uiz^^MiAo}]? (12) 

J a 

By substituting Eq. for u(zq, <po), and making use of 
the identity that 



where 



<% - u(z , M) 



_ e -i*[v-u(za,<i>o)] ^ 



one can then integrate out each A a separately and yield 



In the above 




-iXr/ — \-12 a rnu) a u a (z ,4> )/h 



r, 2 /e 



(14) 



(15) 



For a point r away from the center of the superfluid core, 
along the vector connecting the center to (zo, <fio) on the 
interface, we can calculate the probability W{r) for it to 
be in the polarized fermi gas region: 



/l To 1 
P(r')dr' = - 
-oo 2 



1 + Erf 



r - r 



, (16) 



Therefore the averaged spin polarization is W(r)n p (r) 
and the average pairing gap is (1 — W(r))A s (r), i.e. the 
results of Eq. ([5]). 

Now we consider non-zero temperature. In that case 
we replace Eq. (fl~2|) by 



(17) 



where < ... > denotes the quantum statistical mechani- 
cal average. A convenient way to calculate the average of 
Eq. (|17p is to perform Feynman's path integral in imag- 
inary time. Let A a (r) be the Feynman path of each 
harmonic oscillator in Eq. ([5]), we have 

JYl'VA a (T) 6(n - £ Q X(0KU, M)e~ s 
P(V) = — z (18) 

Here 2, is the partition function, and the action S is given 
by 

S = T f drJ2'[Mrf+^A a (rn (19) 

After using Eq. (fT5|) to represent the delta function, we 
are left with a quadratic path integral which can be per- 
formed exactly to give 



P(v) 



£ 2 7T 



-v 2 /e 



(20) 



~e = — 



ul(z ,a ) 
(ujI + uj 2 a ) ' 



(21) 



and u> n is the bosonic Matsubara frequency. As P(r)) has 
the same form as Eq. (fT4"|) . W(r) will also have the same 
form as Eq. (|12p which leads to the same function form 
of Eq. ©. 

In conclusion, we have demonstrated that even if there 
is no true PP phase in bulk, a PP regime can emerge 
from the quantum and/or thermal fluctuation of the in- 
terface. The width of PP regime depends on the strength 
of the surface fluctuation, and sufficient strong fluctua- 
tion can even lead to disappearance of the equal density 
core, with the pairing gap remaining finite. Though the 
microscopic description of interface fluctuation remains 
to be explored, we have shown that, as long as the theory 
for surface motion is quadratic, there is a universal func- 
tion form for spin polarization and pairing gap, which 
can be compared with existing experiments and tested 
by future experiments as a justification of this scenario. 
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